Anomalous c-axis transport in layered metals 
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Transport in metals with strongly anisotropic single-particle spectrum is studied. Coherent band 
transport in all directions, described by the standard Boltzmann equation, is shown to withstand 
both elastic and inelastic scattering as long as Eft ^> 1. A model of phonon-assisted tunneling via 
resonant states located in between the layers is suggested to explain a non-monotonic temperature 
dependence of the c-axis resistivity observed in experiments. 

PACS numbers: 72.10.-d,72.10.Di 



Electron transport in layered materials exhibits a num- 
ber of unusual properties. The most striking example is 
a qualitatively different behavior of the in-plane (p a t) 
and out-of-plane (p c ) resistivities: whereas the temper- 
ature dependence of p a b is metallic-like, that of p c is ei- 
ther insulating-like or even non-monotonic. At the level 
of non-interacting electrons, layered systems are metals 
with strongly anisotropic Fermi surfaces. A commonly 
used model is free motion along the planes and nearest- 
neighbor hopping between the planes: 



£k = k^/2m a {, + 2 J (1 — cos kj_d) 



(1) 



where k|| and k± are in the in-plane and c-axis com- 
ponents of momentum, respectively, m a }, is the in-plane 
mass, and d is lattice constant in the c-axis direction. For 
the strongly anisotropic case (J <C Ep), the equipoten- 
tial surfaces are "corrugated cylinders" (see Fig.l). 

If the Hamiltonian consists of the band motion with 
spectrum (JT]) and the interaction of electrons with poten- 
tial disorder as well as with inelastic degrees of freedom, 
e.g., phonons, the Boltzmann equation predicts that the 
conductivities are given by 

crf fc = e 2 u{v a v b T t r), a? = 4eVj 2 d 2 (sin 2 (k±d) r tr ), (2) 

where (...) denotes averaging over the Fermi surface and 
over the thermal (Fermi) distribution, v = m a \, /ird is the 
density of states, and t% t is the transport time, resulting 
from all scattering processes (we set U = fee = 1)- If 
T tr decreases with the temperature, both o~ a b and o~ c are 
expected to decrease with T as well. This is not what 
the experiment shows. 

The c-axis puzzle received a lot of attention in con- 
nection to the HTC materials [l[ , and a non- Fermi-liquid 
nature of these materials was suggested to be responsible 
for the anomalous c-axis transport [2j. However, other 
materials, such as graphite TaS 2 [ij, Sr 2 Ru04 or- 
ganic metals @, etc., behave as canonical Fermi liquids 
in all aspects but the c-axis transport. This suggests that 
the origin of the effect is not related to the specific prop- 
erties of HTC compounds but common for all layered 
materials. A large number of models were proposed to 
explain the c-axis puzzle. Despite this variety, most au- 
thors seem to agree on that the coherent band transport 




FIG. 1: Fermi surface corresponding to Eq.([T]) with Fermi 
velocity vectors at two different points. 



in the c-axis direction is destroyed. Although there is no 
agreement as to what replaces the band transport in the 
"incoherent" regime, the most frequently discussed mech- 
anisms include incoherent tunneling between the layers, 



assisted by either out-of-plane impurities 
by coupling to dissipative environment [lc 



by coupl 

II El. 



3 EH, El or 

and polarons 



The message of this Letter is two-fold. First, we ob- 
serve that neither elastic or inelastic (electron-phonon) 
scattering can destroy band transport even in a strongly 
anisotropic metal as long as the familiar parameter EpT 
is large. Nothing happens to the Boltzmann conductivi- 
ties in Eq.© except for becoming very small at high 
temperatures so that other mechanisms, not included in 
Eq-©? dominate transport. This observation is in agree- 
ment with recent experiment [7[ where a coherent fea- 
ture (angle-dependent magnetoresistance) was observed 
in a supposedly incoherent regime. Second, we propose 
phonon-assisted tunneling through resonant impurities as 
the mechanism competing with the band transport. As 
such tunneling provides an additional channel for trans- 
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port, the total conductivity is 



(3) 



where <j res is the resonant-impurity contribution. Be- 
cause <7 ros increases with the temperature, the band chan- 
nel is short-circuited by the resonant one at high enough 
temperatures Accordingly a c goes through a min- 
imum at a certain temperature (and p c = a^ 1 goes 
through a maximum) . We consider phonon-assisted tun- 
neling through a wide band of resonant levels distributed 
uniformly in space. We show that the non-perturbative 
(in the electron-phonon coupling) version of this the- 
ory is in a quantitative agreement with the experiment 
on Sr 2 Ru04 [f|. Due to a similarity between phonon- 
assisted tunneling and other problems, in which inter- 
action leads to the formation of a cloud surrounding the 
electron (such as polaronic effect and zero bias anomaly), 
many ideas put forward ear lier @, M, ED, H H 0, El 
agree with our picture. Nevertheless, we believe that 
only a combination of resonant impurities and electron- 
phonon interaction solves the puzzle of c-axis resistivity 
and provides a microscopic theory for some of the mech- 
anisms considered in prior work. We begin with the dis- 
cussion of the breakdown (or lack of it thereof) of the 
Boltzmann equation. 

One may wonder whether the band transport along 
the c-axis breaks down because the Anderson localization 
transition occurs in the c-direction whereas the in-plane 
transport remains metallic. This does not happen, how- 
ever, because an electron, encountering an obstacle for 
motion along the c-axis, moves quickly to another point 
in the plane, where such an obstacle is absent. More 
formally, it has been shown the Anderson transition oc- 
curs only simultaneously in all directions 0, [l7, 18 1 and 
only if J is exponentially smaller than 1/t. Therefore, 
localization cannot explain the observed behavior. 

Refs.[3, suggested an idea of the "coherent- 
incoherent crossover" . It implies that the coherent band 
motion breaks down if electrons are scattered faster than 
they tunnel between adjacent layers, i.e., if Jt <C 1. Con- 
sequently, the current in the c-direction is carried via in- 
coherent hops between conducting layers. It was noted 
by a number of authors that the assumption about inco- 
herent nature of the transport does not, by itself, explain 
the difference in temperature dependences of o~ a b arid o~ c 
20l 21 1 : due to conservation of the in-plane momentum, 
o~ c is proportional to r both in the coherent and inco- 
herent regimes. Nevertheless, an issue of the "coherent- 
incoherent crossover" poses a fundamentally important 
question: can scattering destroy band transport only in 
some directions, if the spectrum is anisotropic enough 
|22|? We argue here that this is not the case. 

Since we have already ruled out elastic scattering, this 
leaves inelastic one as a potential culprit. We focus on 
the case of the electron-phonon interaction as a source of 
inelastic scattering. For an isotropic metal, the quantum 



kinetic equation is derived from the Keldysh equations of 
motion for the Green's function via the Prange-Kadanoff 
procedure 23] for any strength of the electron-phonon in- 
teraction. In this Letter, we apply the Prange-Kadanoff 
theory to metals with strongly anisotropic Fermi surfaces, 
such as the one in Fig. 1. We show that, exactly as in 
the isotropic case, the Boltzmann equation holds its stan- 
dard form as long as -Ef7"c-p1i 3* 1. Since this form does 
not change between coherent ( JT G _ p h 3> 1) and incoher- 
ent (Jr e - p h *C 1) regimes, it means that the coherent- 
incoherent crossover is, in fact, absent. 

We adopt the standard Frolich Hamiltonian for 
the deformation-potential interaction with longitudinal 
acoustic phonons (u> q = sq) 



H=y £ka k a k 



k 



-X^9 6 q 6 <l + Yl 9 V^k+qGk + ^-q)^) 

q k,q 



Since tunneling matrix elements are much more sensi- 
tive to the increase in the inter-plane distance than the 
clastic moduli, the anisotropy of phonon spectra in lay- 
ered materials, albeit significant, is still weaker than 
the anisotropy of electron spectra (see, e.g., Ref. [241]). 
Therefore, we treat phonons in the isotropic approxima- 
tion, and assume that the magnitude of the Fermi veloc- 
ity is larger than the speed of sound s. 

For a static and uniform electric field, the Keldysh 
component of the electron's Green function satisfies the 
Dyson equation 
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5G K ]_ + [£ K ,®ReG K ]_) 



= -{[E K ,®A] + -[T,®G K } + ) . 



(5) 



Here L — (<9 t + v • Vr + eE • Vk) is the Liouville op- 
erator, A = i(G R — G A ) is the spectral function, T = 
i (T, R — T, A ) , and <g> denotes the convolution in space and 
time. Thanks to the Migdal theorem, the self-energy does 
not depend on electron's dispersion £k = £k — Ep, and 
Eq.© can be integrated over £k- This results in an equa- 
tion 
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i[Re£* S *]_ = 2K K - 1 -[T,g K ) + (6) 



for the "distribution function" 



(7) 



where h = Vk/ |vk| is a local normal to the Fermi surface. 

We consider a linear dc response, when the self-energy 
is needed only at equilibrium. Within the Migdal theory, 
the Matsubara self-energy is given by a single diagram 
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g 2 (q)G(e-oj,-k-q)D(u J ,q). 
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where the dressed phonon propagator 

D- 1 = Do 1 - g 2 Tl 

is expressed through bare one 

D (uJ,q) = - S 2 q 2 / (uj 2 + s 2 q 2 ) 

and polarization operator II which, for Ep > 2 J, is given 
by its 2D form 



n(w,g) = -i/(l-M/ 



V F Q?\ 



We assume that the electron-phonon vertex decays on 
some scale ko shorter than Fermi momentum (ko <C /cf)- 
This assumption allows one to linearize the dispersion 
£k-q ~ £k — Vk • q and simplifies the analysis without 
changing the results qualitatively. As long as J < Ep, 
we have |vk| ~ kpjm a b ~ Vp, where kp is the radius of 
the cylinder in Fig. 1 for J = 0. Despite the fact that 
the electron velocity does have a small component along 
the c-axis, its in-plane component is large (cf. Fig. 1). 
Since it is the magnitude of Vk that controls the Migdal's 
approximation, the problem reduces to the interaction of 
fast 2D electrons with slow 3D phonons. With these 
simplifications, we find 



Im£ R (e,ft) 



i c 



41 -( \k F 

C N 3 

12(1 -C) 2 ^' 



(8a) 
(8b) 



where £ = vg 2 is a dimensionless coupling constant and 
u>d = skp). We see that, despite the strong anisotropy, 
the self-energy remains local, i.e., independent of £k- 

Vertex renormalization leads to two types of correc- 
tions to the self-energy: those that are proportional to 
the Migdal's parameter (s/vp) and those that are pro- 
portional to ms 2 /e. The second type of corrections inval- 
idates the Migdal's theory for temperatures below ms 2 , 



which is about 1 K in a typical metal. For metals with 
anisotropic spectrum the existence of such a scale is po- 
tentially dangerous, since it is not obvious which of the 
masses (light or heavy) defines this scale. We find that 
the in-plane mass (m a b) controls the vertex renormaliza- 
tion for the nearly cylindrical Fermi surface. This shows 
that the Migdal theory for layered metals has the same 
range of applicability as for isotropic metals j25|. 

The rest of the derivation proceeds in the same way as 
for the isotropic case [23}, and the resulting Boltzmann 
equation assumes its standard form. Since no assump- 
tion about the relation between T e _ p h and the dwell time 
(1/J) has been made, the conductivities obtained from 
the Boltzmann equation have the same form regardless 
of whether Jr G _ p h is large or small. In other words, there 
is no coherent-incoherent crossover due to inelastic scat- 
tering in an anisotropic metal [2^ ]. 

The situation changes qualitatively if resonant impu- 
rities are present in between the layers. Electrons that 
tunnel through such impurities are moving with the speed 
controlled by the broadening of a resonant level, i.e., 
much slower than speed of sound. For that reason they 
can not be treated within the formalism outlined above 
and require a separate study. 

To evaluate the resonant-impurity contribution to the 
conductivity, we assume that the impurities are randomly 
distributed in space with density nj mp whereas their en- 
ergy levels uniformly distributed over an interval Eb . The 
tunneling conductance of a bilayer junction is 
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dede'W e 



1 \ 9n' 



(9) 



where W tyt > is a transition probability per unit time and 
n e is the Fermi function. To calculate W^y , we use the re- 



sults of Ref. 3^, 31 1 for the probability of phonon-assisted 
tunneling through a single impurity 



/oo i>oo 
dt ie ltl( - e '- e) / dt2dt 3 e l{t2 - t3Ke - lo) ~ r{t2+t3) (10) 
-00 Jo 

x cxp [ - I 1 " e ~ lH + £ltl ( e ~' it2 ~ X ) I' coth (H) + ( e ~ lt3 + ^ + eltl ^ -!)(!- e ** 3 ) 



where a q = — iAq/ ^/puj q , A is the deformation-potential 
constant, Tl and Tr, are tunneling widths of the resonant 
level, r = Tl + Fr, and eb is the energy of a resonant 
level renormalized by the electron-phonon interaction. In 
the limit of no electron-phonon interaction, Eq. (fl"T))) re- 



produces the well-known Breit-Wigner formula. From 
now on, we consider a wide band of resonant levels: 
Eb ^> T ^> T. Averaging Eq. fTU)) over spatial and en- 



ergy positions of resonant levels, one obtains 
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Ores = Col / de 



1 — coth 



IT) 2T s inh 2 fe) 



/(*) = 



du>- 



( 1 — cos(wt ) ) coth fcp) * s i n ( LJ ^) 



.(11) 



Here cr e i is the conductivity due to elastic resonant tun- 
neling and A = K 2 lo 2 d / ps 5 n 2 is the dimensionless cou- 
pling constant for localized electrons. In the absence 
of electron-phonon interaction, a Ies is temperature inde- 
pendent and given by a e \ ~ ire 2 Tini mp aod/ Eb\<S^, where 
ao is the localization radius of a resonant state and 
rj ~ e e~ d / a ° is its typical width. We note that the 
electron-phonon interaction is much stronger for localized 
electrons than for band ones: X/( ~ (k F d) (vp/s) 3> 1. 
Since typically £ ~ 1, one needs to consider a non- 
perturbative regime of phonon-assisted tunneling. In 
that case, resonant tunneling is exponentially suppressed 
at T = 0: cr res (T = 0) = a ei e- x / 2 . At finite T, we find 



-A/2 



OeZ 



1 + ^ ("^ 

3 V 



1 ^ x/A ' 



As T increases, <7 res growth, resembling the zero-bias 
anomaly in disordered metals and Mossbauer effect. At 
high temperatures (T 3> Au>£>) Ores approaches the non- 
interacting value {a ei). The asymptotic regimes in the 
interval ujjj/^/X <C T <C Acjd can also be studied but we 
will not pause for this here. Notice that, in contrast to 
the phenomenological model of Ref. [|| , there is no simple 
relation between the T-dependences of aj? and cr res . 

To compare our model with the experiment, we extract 
erf from the low-temperature (between 10 and 50 K) o 
axis resistivity of Sr2Ru04 and extrapolate it to higher 
temperatures [j| . The resonant part of the conductivity 
is calculated numerically using Eq. fTTj) . The fit to the 
data for o ei = 43 • 10 3 Vr 1 cm" 1 , uj d = 41 K and A = 16 
is shown in Fig. 2. The agreement between the theory 
and experiment is quite good and the values of the fitting 
parameters are reasonable. An immediate consequence of 
our model is the sample-to-sample variation of the c-axis 
conductivity. Among the layered materials, the largest 
amount of data is collected for graphite 0] • Even within 
the group of samples with comparable in-plane mobili- 
ties, the temperature of the maximum in p c varies from 
40K to 300 K d[33|]. 

To conclude, we have shown that the Boltzmann 
equation and its consequences are no less robust for 
anisotropic metals than they are for isotropic ones. The 
only condition controlling the validity of the Boltzmann 
equation is the large value of Epr, regardless of whether 
r comes from elastic or inelastic scattering. Out-of- 
plane localized states change the c-axis transport rad- 
ically while playing only minor role for the in-plane 
one. While p a b remains metallic, an interplay between 
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FIG. 2: p c vs temperature. Solid: experimental data on 
S^RuCU; dashed: fit into the phonon-assisted tunneling 
model in the non-perturbative regime, Eq. (|ll|l 



phonon-assisted tunneling and conventional momentum 
relaxation causes insulating or non-monotonic depen- 
dence of p c on temperature. This model is in a good 
agreement with the experimental data on S^RuO^ 
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